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ABSTRACT 

Recent  work  developed  a  novel  method  for  determining  tracking  thresholds  for  the  Maximum  Likelihood  Prob¬ 
abilistic  Multi-Hypothesis  Tracker  (ML-PMHT).  Under  certain  "ideal"  conditions,  probability  density  functions 
(PDFs)  for  the  peak  points  in  the  ML-PMHT  log-likelihood  ratio  (LLR)  due  to  just  clutter  measurements  could 
be  calculated.  Analysis  of  these  clutter-induced  peak  PDFs  allowed  for  the  calculation  of  tracking  thresholds, 
which  previously  had  to  be  done  with  time-consuming  Monte  Carlo  simulations.  However,  this  work  was  done 
for  a  very  specific  case:  the  amplitudes  of  both  target  and  clutter  measurements  followed  Rayleigh  distributions. 
The  Rayleigh  distribution  is  a  very  light-tailed  distribution,  and  it  can  be  overly  optimistic  in  predicting  that 
high-SNR  measurements  are  target-originated.  This  work  examines  the  case  where  the  clutter  amplitudes  do 
not  follow  a  Rayleigh  distribution  at  all,  but  instead  follow  a  K-distribution,  which  more  accurately  describes 
active  acoustic  clutter.  This  will  provide  a  framework  for  determining  accurate  tracking  thresholds  for  the 
ML-PMHT  algorithm. 

Keywords:  ML-PMHT,  Multistatic,  Bistatic,  Active  Sonar,  Tracking,  Threshold  Determination,  K-distribution 

1.  INTRODUCTION  AND  ML-PMHT  BACKGROUND 

The  Maximum  Likelihood  Probabilistic  Multi-Hypothesis  Tracker  (ML-PMHT)  is  a  powerful  tracking  algorithm 
that  can  be  implemented  in  a  multistatic  active  sonar  framework.  It  has  its  origin  in  a  similar  algorithm,  the 
Maximum  Likelihood  Probabilistic  Data  Association  (ML-PDA)  tracker,  which  was  first  developed  in  [16]  and 
was  expanded  in  [17]  and  [9].  The  ML-PMHT  algorithm  followed  from  these,  as  well  as  the  work  in  [5],  [23], 
[24],  and  [25].  The  development  of  the  ML-PMHT  algorithm  as  a  multistatic  sonar  algorithm  is  described  in 
detail  in  [22]  —  it  is  only  briefly  summarized  here. 

At  its  heart,  ML-PMHT  is  a  very  simple  algorithm.  Several  assumptions  about  the  target  and  the  environ¬ 
ment  are  made  [6],  and  using  these  assumptions,  a  log-likelihood  ratio  is  constructed.  This  log-likelihood  ratio 
for  a  single  target  is  given  by 


Nw  rrii  s  >. 

A(x>^)  =5/X]ln{l  +  ^kpj(i)p[zi(f)|x]|  (1) 

i=  1  j— 1  ^  ' 

where  x  is  a  target  parameter  vector,  Z  is  the  entire  set  of  measurements  in  a  batch  of  scans,  Nw  is  the  number 
of  scans  in  the  batch,  to,;  is  the  number  of  measurements  in  the  Ith  scan,  7Ti  is  the  prior  probability  that  any 
given  measurement  is  from  the  target,  ttq  is  the  prior  probability  that  any  given  measurement  is  from  clutter, 
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V  is  the  measurement  search  volume,  pj  (i)  is  the  amplitude  likelihood  ratio  for  the  jth  measurement  in  the  ith 
scan,  z j(i)  is  the  jth  measurement  in  the  ith  scan,  and  p[zj(i)\x\  is  a  target-centered  Gaussian.  It  is  stressed  that 
ML-PMHT  is  a  batch  tracker,  as  opposed  to  a  recursive  tracker  such  a  Kalman  filter,  and  it  is  non-Bayesian  — 
it  assumes  that  the  target  follows  some  deterministic  (unknown)  parameterizable  trajectory  which  is  corrupted 
by  measurement  noise. 

To  apply  ML-PMHT  as  a  multistatic  tracker,  a  batch  of  measurements  is  accumulated  and  inserted  into  (1), 
and  the  LLR  is  optimized  over  the  target  parameter  vector  x.  If  the  maximum  LLR  value  found  is  above  a 
certain  threshold,  a  target  is  declared;  if  the  maximum  LLR  value  is  below  this  threshold,  a  no-target  condition 
is  assumed.  A  problematic  part  of  actually  implementing  ML-PMHT  is  determining  exactly  what  this  thresh¬ 
old  should  be.  In  some  ML-PDA/ ML-PMHT  implementations  [14],  [19],  and  [20],  “optimal"  thresholds  were 
determined  mainly  through  trial  and  error  (which  obviously  predicates  use  of  the  tracker  in  real-time).  Work 
in  [7]  did  develop  a  method  for  determining  thresholds  empirically  via  Monte  Carlo  testing,  but  this  approach 
is  very  computationally  intensive  and  time  consuming;  again,  it  is  not  possible  to  use  this  method  of  threshold 
determination  with  real-time  processing.  A  fast  and  accurate  method  for  ML-PMHT  threshold  determination 
is  required. 


2.  DETERMINING  ML-PMHT  TRACKING  THRESHOLDS 

Recent  work  in  [21]  developed  a  method  for  determining  ML-PMHT  tracking  thresholds  theoretically  (and 
rapidly),  but  it  was  limited  to  only  several  specific  types  of  clutter  measurements.  This  work  is  summarized 
below  and  is  then  expanded  to  a  more  general  framework  that  numerically  computes  ML-PMHT  tracking 
thresholds  for  a  wider  variety  of  cases.  The  resultant  overall  framework  for  determining  clutter  thresholds  is 
shown  in  Figure  1. 

2.1  Determining  peak  ML-PMHT  LLR  PDFs 

The  clutter  threshold  determination  starts  with  an  individual  clutter  measurement.  In  this  work,  the  focus  is  on 
3-dimensional  measurements.  The  individual  measurements  are  azimuth,  time-delay,  and  range  rate.  Consider 
(1)  for  just  a  single  measurement  This  can  be  written  as 


Ai(x,  z)  =  In 


1  +  Ke 


Pi —mi y 


I  P2-M2) 

+ - 9^ 


I  U3-M3)' 


(2) 


where 

™  (3) 

7T0  V87T3  (71(720-3 

Here,  pi,  /i->  and  //,3  are  the  mean  target  locations  in  measurement  space,  a\,  cr2  and  03  are  the  respective  mea¬ 
surement  errors,  and  Vi  ,  V2,  and  V3  are  the  individual-dimension  measurement  search  volumes.  For  the  time 
being,  it  is  assumed  that  no  amplitude  information  is  being  processed,  so  pj(i)  =  1. 

One  of  the  assumptions  of  the  ML-PMHT  algorithm  is  that  clutter  measurements  are  uniformly  distributed 
in  measurement  space;  with  this,  (2)  can  be  treated  as  simply  a  random  variable  (RV)  transformation;  i.e.  if  Z\, 
Z2,  and  Z’i  are  uniformly  distributed  RVs  in  azimuthal,  time-delay,  and  range-rate  measurement  space,  then  W 
is  simply  a  new  transformed  RV,  where 

W  =  Ai(x,  z)  (4) 

Using  standard  random  variable  transformation  techniques,  then,  it  is  possible  to  obtain  the  PDF  for  W  under 
various  cases.  This  "single-measurement  PDF"  of  W  represents  the  range  of  values  (with  their  correspond¬ 
ing  probability  densities)  that  the  ML-PMHT  LLR  can  take  on  from  a  single  clutter  measurement.  In  [21], 
analytic  expressions  for  the  single-measurement  PDFs  were  derived  for  the  following  cases:  1-dimensional 
measurements  (time-delay  only  or  azimuth  only),  2-dimensional  measurements  (azimuth  and  time-delay), 
3-dimensional  measurements  (azimuth,  time-delay,  and  range-rate),  and  2-dimensional  measurements  with 
Rayleigh-distributed  amplitudes  for  both  target-  and  clutter-originated  measurements. 
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Figure  1.  Framework  for  determining  clutter  peak  PDF 


ML-PMHT  is  a  batch  tracker;  it  is  necessary  to  obtain  the  PDF  of  a  batch  of  measurements,  instead  of  just 
the  PDF  of  a  single  measurement.  Fortunately,  another  one  of  the  ML-PMHT  assumptions  (and  this  can  be  seen 
in  the  structure  of  the  ML-PMHT  LLR)  is  that  measurements  at  different  times  are  assumed  to  be  independent. 
Thus,  the  single-measurement  RVs  resulting  from  the  transformation  described  by  (4)  will  also  be  independent 
(and  identically  distributed).  As  a  result,  the  PDF  of  a  batch  of  N  measurements  can  be  computed  simply  by 
convolving  the  single-measurement  PDF  of  W  N  —  1  times  (the  PDF  of  a  sum  of  two  independent  RVs  is  just 
the  convolution  of  the  individual  RV  PDFs).  This  "batch  PDF"  probabilistically  describes  the  range  of  possible 
values  of  the  ML-PMHT  LLR  (1). 

This  batch  PDF  describes  all  possible  values  for  the  ML-PMHT  LLR;  we  are  concerned  with  the  peak  value 
of  the  ML-PMHT  LLR.  To  obtain  the  PDF  for  this  peak  value,  extreme-value  (EV)  theory  is  used  [8],  [10],  [11], 
[12],  [13],  [15].  According  to  EV  theory,  in  this  case,  the  PDF  of  the  peak  value  of  the  LLR  can  be  approximated 
with  a  Gumbel  distribution,  which  is  expressed  as 


f(x)  =  p  exP 


(5) 


where  v  is  the  location  parameter  and  /?  is  the  shape  parameter.  In  an  "actual"  ML-PMHT  implementation, 
a  non-linear  numerical  optimization  scheme  would  be  run  to  optimize  (1)  over  x.  In  theory,  if  the  batch  ML- 
PMHT  LLR  is  sampled  M  times,  where  M  is  some  large  number,  the  value  of  the  maximum  sample  should 
be  the  same  as  the  maximum  LLR  value  determined  by  numerical  optimization  (i.e.  both  results  would  come 
equally  close  to  the  "true"  LLR  maximum).  With  that  idea  in  mind,  if  the  accuracy  of  the  actual  optimizer  is 
known,  then  is  possible  to  compute  a  value  for  M  that  would  produce  (via  sampling)  a  maximum  with  this 
same  accuracy.  From  extreme-value  theory,  the  PDF  of  the  maximum  sample  from  a  series  of  length  M  IID  RV 
samples  takes  on  a  Gumbel  distribution.  The  parameters  v  and  /3  in  (5)  are  a  function  of  M  and  the  underlying 
cumulative  distribution  function  (CDF)  of  the  sampled  RV.  The  key  point  is  that  the  numerical  optimization  is 
represented  as  taking  the  maximum  from  a  large  sample  of  IID  RVs;  we  compute  how  many  times  it  is  necessary 
to  sample  the  batch  PDF  to  obtain  the  same  accuracy  as  the  numerical  optimization  achieves  (note  that  actual 
optimization  via  sampling  is  not  done),  and  then  this  value  M  is  used  along  with  the  CDF  of  the  underlying 
batch  RV  to  determine  the  "peak"  PDF  —  this  statistically  characterizes  the  peak  point  in  the  ML-PMHT  LLR. 
Finally,  the  Neyman-Pearson  Lemma  from  classical  detection  theory  [18]  is  used  —  a  desired  probability  of  false 
track  ( Pft )  is  selected,  which  in  turn  determines  the  tracking  threshold. 

One  of  the  steps  in  this  process  is  to  perform  the  transformation  of  a  uniformly  distributed  RV  to  the  single¬ 
measurement  RV  using  the  single-measurement  ML-PMHT  LLR  as  the  transforming  function.  Again,  this 
was  done  analytically  in  [21]  for  four  cases:  1-dimensional  measurements,  2-dimensional  measurements,  3- 
dimensional  measurements,  and  2-dimensional  measurements  with  a  Rayleigh  amplitude  ratio.  It  does  not 
appear  to  be  possible  to  analytically  calculate  the  expression  for  a  single-measurement  transformed  RV  for  the 
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2-dimensional  case  with  a  clutter  amplitude  distribution  other  than  Rayleigh,  nor  does  it  appear  to  be  possible 
to  calculate  the  transformed  PDF  for  any  amplitude  distribution  in  the  3-dimensional  case  (the  integrals  for  the 
probability  transformations  are  intractable  analytically).1  In  this  work,  numerical  methods  are  used  to  compute 
the  single-measurement  PDFs  for  the  cases  where  analytic  results  are  not  possible.  We  are  particularly  interested 
in  cases  where  the  clutter  is  modeled  as  being  K-distributed,  which  is  a  more  accurate  model  of  "real"  clutter 
[3]. 

2.2  Amplitude  likelihood  ratios 

Previous  work  dealing  with  amplitudes  [17]  assumed  that  both  the  clutter  and  the  target  amplitude  followed 
a  Rayleigh  distribution.  In  intensity  space  (amplitude  squared)  these  become  exponential  distributions.  (From 
here  on,  work  is  done  in  intensity  space  because  the  exponential  distribution  is  very  simple  to  work  with.) 
Using  this  model,  clutter  intensity  is  modeled  as  a  unit  mean  exponential,  and  the  target  intensity  is  modeled  as 
an  exponential  with  mean  a2.  Additionally,  both  the  target  and  clutter  distributions  are  thresholded  at  certain 
level  r,  which  is  the  sonar  detector  threshold.  The  amplitude  LR  that  results  is 


Pj(i)  =  f{v) 


1  e 

a2  eT~v 


V  >  T 


(6) 


where  v  is  the  intensity,  a2  is  the  expected  intensity,  and  r  is  the  detector  threshold,  in  units  of  intensity.  (It  is 
important  to  note  that  regardless  of  whether  measurements  are  treated  in  amplitude  space  or  intensity  space, 
the  likelihood  ratios  come  out  the  same.) 

A  "problem"  with  the  Rayleigh/ exponential  distributions  when  they  are  used  to  describe  clutter  is  that  they 
are  very  "light-tailed"  distributions.  The  result  of  this  light  tail  is  that  the  amplitude  LLR  for  high-SNR  mea¬ 
surements  becomes  very  large  —  this  is  shown  in  Figure  3.  When  a  large-amplitude  measurement  is  plugged 
into  pj(i)  in  (1),  it  essentially  dominates  all  other  measurements.  As  a  result,  the  ML-PMHT  LLR  gets  all  of  its 
value  (and  the  resulting  discovered  target  gets  its  positional  measurement)  from  this  single  measurement.  This 
is  bad  enough  if  the  measurement  is  actually  from  the  target  (it  is  difficult  to  get  an  accurate  target  solution  off 
of  one  measurement);  it  is  even  worse  if  the  high-amplitude  measurement  is  actually  from  clutter.  In  such  a 
case,  a  false  track  is  practically  guaranteed. 

Recent  work  [1],  [2],  and  [3]  has  shown  that  the  K-distribution  provides  a  more  accurate  model  of  active 
acoustic  clutter.  The  probability  distribution  for  the  clutter  intensity  in  this  case  is  given  by 


/  0) 


2 

ATK) 


{ak- 1)/2 

Kuk-i 


v  >  0 


(7) 


Flere,  ak  is  a  K-distribution  parameter,  A k  -  1  /oik,  and  K„  is  the  Basset  function  (a  modified  Bessel  function  of 
the  second  kind)  [4].  The  parameter  ak  for  this  distribution  determines  the  "heaviness"  of  the  tail;  the  smaller 
the  value  of  ak,  the  heavier  the  tail  (i.e.  the  more  the  probability  mass  is  spread  out  to  the  right).  This  is  shown 
in  Figure  2  —  for  large  values  of  ak,  the  distribution  approaches  that  of  the  unit  exponential;  while  for  small 
values  of  ak,  the  tails  of  the  distributions  become  much  heaver. 

Such  a  distribution  has  a  very  interesting  effect  on  the  amplitude/intensity  LLR  (it  is  important  to  keep  in 
mind  for  the  following  discussion  that  the  target  amplitude/ intensity  is  still  modeled  as  a  Rayleigh/ exponential 
RV  with  parameter  a2).  The  amplitude  LLR  using  K-distributed  clutter  is  shown  in  Figure  4.  For  small  SNR 
values,  the  LLR  is  negative,  which  indicates  that  the  measurement  is  more  likely  to  have  been  clutter-originated. 
For  medium  SNR  values  (10-23  dB  in  this  case),  the  LLR  is  positive,  indicating  that  the  measurement  is  more 
likely  to  have  been  target-originated.  However,  for  large  SNR  measurements,  the  LLR  is  again  negative,  which 
says  that  the  measurement  is  most  likely  clutter-originated.  This  stands  in  stark  contrast  to  the  Rayleigh  model, 
which  mathematically  says  that  large  SNR  measurements  come  from  the  target,  essentially  without  a  doubt. 

1It  was  not  proven  that  these  cases  do  not  have  analytic  solutions.  However,  the  authors  were  unable  after  much  effort  to 
analytically  solve  the  integrals  necessary  for  the  probability  transformations. 
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SNR  (dB) 

Figure  2.  Squared  K-distribution  PDF  and  unit  exponential  PDF 


SNR  (dB) 

Figure  3.  LLR  for  Rayleigh  distribution  clutter 


Figure  4.  LLR  for  K-distribution  clutter 


2.3  Numerically  computing  single-measurement  LLRs 

At  this  point,  it  is  necessary  to  develop  a  framework  for  numerically  computing  the  single-measurement  PDFs 
that  result  from  transforming  a  uniform  random  variable  by  the  single-measurement  ML-PMHT  LLR  (2).  This 
is  done  by  rewriting  the  single-measurement  ML-PMHT  LLR  in  a  "random  variable"  form 


W  =  ln(l  +  Z) 


(8) 


where 


Z  =  XY 


(9) 


Here  W ,  Z,  X,  and  Y  are  all  random  variables;  we  are  ultimately  trying  to  compute  the  PDF  of  W.  Writing  the 
expression  for  Z  gives 


Z  =  pj(i)  Ke 
x 


_  \  (zj-AM)  |  (>2-/-*2)  |  (z^-lizY 


(10) 


It  is  actually  possible  to  get  an  analytic  solution  for  the  PDF  of  Y  using  the  standard  one-to-one  RV  transforma¬ 
tion  technique.  Start  by  defining  the  exponential  term  as  single,  transformed  random  variable  —  a  sum  of  three 
shifted,  scaled,  and  squared  uniform  random  variables 


U  = 


{zi  -  Mi)2  .  (Z2  -  M2)2  .  (-3  -  M3)2 


2  tj\ 


2ct22 


The  PDF  for  U  comes  out  to  be2 


2  °i 


(11) 


/  \  o  °T  <72  <73 
Pu(U)  =  27T-— — VM 


V1V2V3 

This  in  turn  leads  to  an  expression  for  the  PDF  of  Y 


0  <  u  <  min( 


(12) 


,  .  0  cricr2cr3  /  A 

pY{y)  =  2tt  2  In  - 

V1V2V3  V  y  y 


0  <  y  <  K 


(13) 


The  PDF  for  X  cannot  be  determined  analytically.  To  calculate  it  numerically,  start  with  a  squared  K- 
distributed  random  variable  with  parameter  ctk,  and  transform  it  with  the  amplitude/intensity  LLR 


f(v)  = 


dy  exp[(r  —  v)/cr2} 


CJrf^T)  (  a”  ) 


(afc-l)/2 


K, 


OCk  —  1 


/  9  /  v  \ 

V  V ) 


(14) 


Here,  C  is  a  numerically-calculated  normalization  constant  for  a  thresholded  squared  K-distribution  PDF. 
Again,  the  transformation  is  performed  with  the  standard  one-to-one  RV  transformation  method.  Results  for 
the  PDF  of  X  are  shown  in  Figures  5,  6,  and  7  for  three  different  values  of  au-  (Because  some  of  the  theoretical 
PDFs  for  X  looked  very  "strange"  —  see  Figure  5  as  an  example  —  all  of  these  PDFs  were  verified  correct  with 
empirical  testing.) 

With  the  PDFs  of  X  and  Y  now  available,  the  PDF  of  Z  —  XY  can  be  computed  by  introducing  the  aux¬ 
iliary  random  variable  U  —  Z ,  performing  a  standard  two-to-two  RV  transformation  to  get  the  bivariate  PDF 
Pz,u(zi  u)>  and  then  integrating  out  U  to  get  Pz(z).  Finally,  the  PDF  of  W  is  obtained  with  one  final  one-to-one 
numerical  transform  described  by  (4).  Resultant  PDFs  for  W  for  three  values  of  a k  are  shown  together  in  Figure 
8.  (Again,  note  that  these  and  all  following  results  are  for  3-dimensional  measurements.) 


Probability  density 
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Figure  5.  PDF  of  X  for  au  =  0.1  Figure  6.  PDF  of  X  for  ctk  =  0.5  Figure  7.  PDF  of  X  for  oik  =  5 


X  10“4 


Figure  8.  Single-measurement  W  PDFs  for  ctk  =  (0.1,  0.5,  5.0} 


Figure  9.  Batch  PDFs,  peak  PDFs,  Figure  10.  Batch  PDFs,  peak  PDFs,  Figure  11.  Batch  PDFs,  peak  PDFs, 

and  clutter  thresholds  for  a  =  0.1  and  clutter  thresholds  for  a  =  0.5  and  clutter  thresholds  for  a  =  5 
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Table  1.  Comparison  of  model  vs.  empirical  tracking  thresholds  (dimensionless  units  for  LLR).  Model  thresholds  are  from 
the  method  developed  in  this  work,  while  empirical  thresholds  and  confidence  intervals  are  from  Monte  Carlo  testing. 


Olk 

Model 

Empirical 

Empirical  95% 
Confidence  Interval 

0.1 

40.2 

39.5 

[38.7, 40.3] 

0.5 

38.1 

37.5 

[37.0, 38.1] 

0.8 

32.2 

32.9 

[32.4, 33.4] 

1.0 

30.7 

31.2 

[30.5, 32.1] 

2.0 

28.2 

28.2 

[27.5, 29.0] 

5.0 

24.5 

24.3 

[23.7, 24.8] 

10 

22.1 

22.8 

[22.2, 23.3] 

50 

19.0 

20.5 

[20.1,20.8] 

Rayleigh 

18.5 

20.0 

[19.8, 20.3] 

3.  RESULTS 

With  the  single-measurement  PDFs  calculated  numerically,  the  remainder  of  the  framework  shown  in  Figure 
1  is  applied.  Resultant  batch  PDFs,  extreme-value  (peak)  PDFs  and  thresholds  (for  Pft  =  0.01)  are  shown 
in  Figures  9,  10,  and  11.  All  of  the  calculated  model  thresholds  are  then  compared  to  empirically  determined 
thresholds,  and  in  almost  all  cases  the  results  agree  extremely  well  with  the  empirical  values  —  in  most  cases 
the  model  values  fall  within  the  95  percent  empirical  error  bounds,  and  for  instances  where  this  is  not  the  case 
the  values  still  are  reasonably  close.  Also  note  that  the  results  are  consistent  with  what  is  expected  —  as  oik  gets 
larger,  the  tail  of  the  K-distribution  gets  lighter  and  lighter,  approaching  that  of  a  Rayleigh.  This  exactly  what  is 
seen  in  Table  1.  The  results  for  higher  values  of  a.k  approach  those  of  Rayleigh-distributed  clutter. 

The  empirical  values  in  Table  1  are  "correct."  They  are  the  result  of  many  clutter-only  simulations  that  are 
run  through  the  ML-PMHT  tracker  and  then  optimized  in  a  Monte-Carlo  fashion,  but  they  are  also  compu¬ 
tationally  intensive  to  calculate.  Each  of  the  empirical  thresholds  in  Table  1  took  several  hours  of  run-time 
to  calculate.  In  contrast,  the  model  values  each  take  only  several  seconds  to  compute,  and  they  are  accurate 
enough  to  use  in  the  ML-PMHT  tracking  framework. 

4.  CONCLUSION 

We  have  expanded  on  the  methods  developed  in  [21]  to  develop  ML-PMHT  thresholds  for  K-distributed  clutter, 
which  is  a  more  realistic  distribution  for  describing  active  clutter  amplitude/intensity  measurements.  Trans¬ 
formed  single-measurement  PDFs  were  numerically  calculated  for  various  values  of  ak,  the  parameter  that  con¬ 
trols  the  K-distribution  shape.  These  numerical  PDFs  were  then  processed  through  the  previously  developed 
framework  to  obtain  batch  PDFs,  peak  PDFs  and  finally  ML-PMHT  clutter  thresholds.  The  ability  to  rapidly 
and  accurately  determine  these  clutter  thresholds  for  a  more  realistic  clutter  amplitude  parameterization  will 
improve  ML-PMHT  as  a  multistatic  active  tracker. 
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